Rabi spectroscopy for a two-level system 
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We have analyzed the interaction of a dissipative two level quantum system with high and low 
frequency excitation. The system is continuously and simultaneously irradiated by these two waves. 
If the frequency of the first signal is close to the level separation and the second one is tuned to the 
Rabi frequency, it is shown that the response of the system exhibits an undamped low frequency 
oscillation. The method can be useful for low frequency Rabi spectroscopy in various physical 
systems which are described by a two level Hamiltonian, such as nuclei spins in NMR, double well 
quantum dots, superconducting flux and charge qubits, etc. As an example, the application of the 
method to the readout of a flux qubit is considered. 
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It is well known that under resonant irradiation, a 
quantum two level system (TLS) can undergo coherent 
(Rabi) oscillations. Such oscillations occur if microwaves 
are in resonance with the spacing between the energy 
levels. In this case the level occupation probability will 
oscillate with a frequency proportional to the amplitude 
of resonant field [l| . The effect is widely used in molecu- 
lar beam spectroscopy 2] , and in quantum optics |3( • 

Normally, such oscillations are damped out with a rate, 
which is dependent on how strongly the system is cou- 
pled to the environment. In recent years, progress in the 
investigations of scalable solid state qubits has necessi- 
tated the development of different methods of advanced 
qubit control. One of the proposals in this direction was 
to maintain the Rabi oscillations for arbitrary long times 
by the synchronization of the phase of Rabi oscillations, 
with that of the ideal solid state detector with the aid of 
a quantum feedback loop 

In this work we propose another method for maintain- 
ing such oscillations. We consider a TLS which is irra- 
diated continuously by two external sources. The first 
one with a frequency uq, which is close to the energy gap 
between the two levels, excites the low frequency Rabi 
oscillations. A second, low frequency source tuned to the 
Rabi frequency causes the oscillations to persist. 

We start with a Hamiltonian of a TLS subjected to a 
both high and low frequency excitation: 
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a z F cos uj t — a z G(t). (1) 



Hamiltonian JJ) is written in the localized state basis. 
Here A is the tunnelling rate between localized states, 
e is the bias, F is the coupling between the TLS and 
high frequency field, and G(t) is a low frequency exter- 
nal force. In the eigenstates basis, which we denote by 
upper case subscripts for the Pauli matrices ax , ay , <?z > 



Hamiltonian 0} reads: 
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where A e = y/ A 2 + £ 2 is the gap between two energy 
states. 

A more realistic description of the TLS requires the 
inclusion of the dissipative environment in Hamiltonian 
JJJ, which for a spin-boson model of the bath results 
in the Bloch-Redfield equations for the matrix operators 
ax, 0y, az, @, 0] • For weak driving (/ -C A), and for 
weak coupling of the TLS to the bath these equations 
can be approximated by Bloch-type equations 



(a z ) = (2 f cos io t + 2g(t))(a Y ) -T z ((a z ) - Z ) , (3) 



(a Y ) = - (2/ cos u t + 2g(t)) (a z ) 

(a X ) -T (ay), (4) 
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where / = AF/fiA e , g(t) = AG(t)/?zA £ , and Z = 
— tanh (A e /fcsT) is the equilibrium polarization of the 
system in the absence of external excitation sources 
(/ = 0, g = 0). The angled brackets in Eqs. ©, QJ, 
and JSJ denote the trace over reduced density matrix 
p(t), which is obtained by tracing out all environment 
degrees of freedom: (ax) — Tr(axp(t)), etc. In order to 
simplify the problem we assume the relaxation, T z and 
dephasing, T, rates in ar e time- independent, i.e. 

the rates are slowly varying functions on the scale of Rabi 
period which is on the order of 1//. 
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The desired solution of Eqs. ©, ®, and JSJ is: 

(<rz) = Z(t), 



(6) 



(cry) = Y(t) + A(t) cos(w i) + B{t) sin(w t), (7) 



(<jx) = X(t) + C(t) cos{uj t) + D(t) sm(wo*), (8) 

where the quantities X, Y, Z, A, B, C, and D are slowly 
varying, compared to the high frequency ujq operators. In 
the rotating wave approximation we obtain from Eqs. Jjyl, 
(@J , and the response of the system to a small external 



force g(t). The Fourier components for the slowly varying 
A(ui), and B{uj), and for the low frequency persistent 
response of the TLS to a small low frequency excitation, 
Z(uj), Y(u>), and X(u>) are as follows: C(u) = B(u>), 
D(co) = -A(w), 



Z w) = g(uj)—f PzZq — -— , 
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where 6 is a high frequency detuning, 5 = ujq — A £ /h, £Ir 
is the Rabi frequency, CIr = \/ S 2 + f 2 , s(u>) — (CIr — u) + 
iT)(Q R + uj- %T)(iu) + T z ) - f 2 (T z - T), and 



Tz(T 2 + 5 2 ) + f 2 V v ; 

PzZq is the noncquilibrium polarization, i. e., the steady 
state population difference between the qubit energy lev- 
els in the case when the high frequency excitation is ap- 
plied to the qubit. When deriving the expressions ©- 
(|13J) we assume the force g to be small, thereby we ne- 
glect the higher order terms of g. In addition, we keep 
only the terms which oscillate within the bandwidth of 
the Rabi frequency and neglect the terms which are of 
the order of ft//A ej ar/A £ , hT z /A £ . 

As is seen from these expressions the persistent low 
frequency oscillations of the spin components Z(i), X(t), 
Y(t) appear as the response to the low frequency external 
force, git), only in the presence of the high frequency 
excitation (/ ^ 0). 

Exactly at resonance (S = 0) A(uS) — 0, Z(u>) — 0, but 
B(u>) ^ 0. Hence, as is seen from Eqs. ©, 0, ®, in 
this case the population of the two levels are equalized, 
spin circularly rotates in the XY plane with frequency 



luq, with the center of the circle being precessed with the 
Rabi frequency. 

The key point of the method we described above is 
that it allows for the detection of the high frequency re- 
sponse of a TLS at a frequency which is much less than 
the gap frequency. The low frequency dynamics of the 
quantities (<Jx) , (cy), (<?z) bears information about the 
relaxation, Tz, and dephasing, T, rates. The experimen- 
tal realization of this method requires the measurement 
of the difference of low frequency response of the system 
of interest with and without high frequency excitation. 

The detection scheme depends on the problem under 
investigation. For example, the method can easily be 
adapted for NMR. As is seen from 0J, the NMR case 
corresponds to the polarization of the sample with the 
fields A and e along the x and z axes, respectively, with a 
high frequency excitation Fcos(u>ot) and a low frequency 
probe G(t) being applied along the z axis. 

Another example are the macroscopic quantun TLS, 
such as superconducting qubits based on mesoscopic 
Josephson junctions Q For these qubits the Rabi os- 
cillations have been detected through the statistics of 
switching events 0, IToL fill H^ ] . In parallel, the detec- 
tion of coherent oscillations through a weak continuous 
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measurement which does not destroy the quantum coher- 
ence of the TLS was proposed [lj, [lj, [l5j . Generally, the 
signature of such oscillations can be found in the noise 
spectrum of the TLS, which was recently demonstrated 
for a superconducting flux qubit 0]. Let us consider 
last example in more detail, namely, the application of 
the method to a persistent current qubit, which is a su- 
perconducting loop interrupted by three Josephson junc- 
tions 0], 01- The basis qubit states have an opposing 
persistent current and the operator of the persistent cur- 
rent in the qubit loop reads I q = I q a z . In eigenstate basis 
the average current, (I q ) is: 



(15) 



This current can be detected through the variation of 
its magnetic flux either by a DC SQUID 0] or by a 
high quality resonant tank circuit inductively coupled to 
the qubit [iJ IU . 

For the flux qubit the bias e is controlled by a dc flux 
$x: e = Ejfx, where fx = ($x/$o - 1/2), Ej = 
§oI q /2ir, and <I>o = h/2e is the flux quantum. The qubit 
is inductively coupled to a high quality resonant tank 
circuit with inductance Lt, capacitance Ct, and quality 
factor Qt- The mutual inductance between the qubit and 
the tank M — k{L q LT) x ^ 2 , where k is the dimensionless 



coupling parameter, and L q is the inductance of the qubit 
loop. The tank is biased by a low (compared to the 
gap) frequency current It, = Iocos(ujt). This readout 
circuit has proven to be successful for the investigation 
of quantum properties of the flux qubit 0, 0, 12^ . 
The qubit+tank system is described by the Hamiltonian 
H + Ht, where Ht is the Hamiltonian of the tank, and 
H is given in (JJJ, with G(t) = MI q lT(t), where I q is the 
current flowing in the qubit loop, and It is the current 
in the tank coil. The voltage across the tank, V(t) = 
Vrcos(u)t + x), is given by 20 : 



V + j T V + lo t V = Mw T 



dt 



- LO T 
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(16) 



where jt — <-^t/Qt, and lot — (LtCt)^ 1 ^ 2 is the tank 
resonance frequency, which is tuned to the Rabi fre- 
quency {lot — £Ir). 

From Eqs. |@J|, and J5J we obtain: 



( -eTz [(crz) - Z ] + AT (a x ) 



d(I q 
dt 

In order to get the tank response we have to keep, in the 
right hand side of (jT7fl . only the low frequency parts of 
the averaged Pauli operators: 
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V{io) (lj t - uj 2 + iujr) = M^Il (- e T z Z{io) + ATX(lo) + yAF(u) ) - L / ■/,,. 



(18) 



By taking into account that G(ui) — Mlqlrfa), where It{lo) = —iV(u))/loLt, we obtain with the aid of the expressions 
(JSJ- ljT3)l the low frequency detuning £, and the damping, Ft of the tank: 
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The quantities £ and Tt are related to the voltage amplitude Vt, and the phase, x : Vr = lolo t LtIq/ \/i 2 + ^T 2 -,, 
and tan% = ^/loTt- The functions /i, /2, and d, which account for the resonant properties of the Rabi oscillations, 
are as follows: 



= [(2rr z + lo 2 ) (n% + r 2 - lo 2 ) + 2co 2 r (r z - 2r) - 2/ 2 r (r z - r)] , 



(21) 



/ 2 (w) = [r z (n| + r 2 - co 2 ) + / 2 (r - v z ) - 2r (n* + r 2 + 2tt z )} , 



(22) 



d{u) = 



(n 2 +r 2 -Lo 2 ) 2 + 4uj 2 T 2 



r|) + / 4 (r z - r) 2 - 2/ 2 (r z - r) [r z (n 2 R + r 2 - ^ 2 ) - 2lo 2 t] . (23) 



The equations (|19II20|I are obtained for the case AS/h ^> f 2 . In the opposite case, in particular, for zero high 
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frequency detuning (5 = 0) the result is: 
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where p = p z (s = o) = r z r/(r z r + f). 

In the absence of high frequency excitation (/ = 0) 
we obtain from Eqs. {fU |2UJ| (or from Eqs. 113 a 
purely inductive response (r^ = It, Pz = 1), which de- 
scribes the adiabatic evolution of the qubit in the ground 
state 21]. The experimental study of the adiabatic evo- 
lution provides us with information about the energy gap 
between the two levels 22]. Therefore, the difference 
of the low frequency response of the persistent current 
qubit with and without high frequency excitation can 
provide additional information about the qubit's damp- 
ing rates, T and Tz (see Eqs. EH 12011 . As an illus- 
tration of the method we show in Fig. ^ how the de- 
pendence of the tank phase \ on the flux changes with 
the application of the high frequency excitation. The 
graphs have been calculated for zero high frequency de- 
tuning (5 = 0) from Eq. I|24|) for the following parame- 
ters: u/2ir = u t /2tt = 6 MHz, k = 0.03, Q T = 2000; 
Ej = 1.32 x 10~ 22 J, E c = 2.14 x 10~ 24 J, L q = 40 pH, 
/, = 280 nA, A/2tt = 1 GHz, T z /2ir = 0.1 MHz, 
T/2tt = 4 MHz, and T = 10 mK. At a relatively low 
power of the irradiation the form of the curves remains 
unchanged, with the amplitudes of the dips being condi- 
tioned by the factor P Q in Eq. (2U (Fig. Hk). At higher 
powers the second term in the curled brackets of Eq. (|24H , 
which describes the influence of Rabi oscillations, comes 
into play. As a result, the forms of the phase curves are 
changed drastically (Fig. HJ)). 

In conclusion, we propose a method to maintain the 
Rabi oscillations in a dissipative TLS by irradiating it 
with the two external sources. The first, high frequency 
source, which is tuned to the energy gap of the system, 
excites low frequency Rabi oscillations. The second, low 
frequency source, which is tuned to the frequency of the 
Rabi oscillations, causes the oscillations to persist. When 
applied to a readout of the flux qubit inductively coupled 
to a tank circuit, the method allows for the experimental 
determination of the relaxation and dephasing rates of 
the qubit. 

We are grateful to D. Averin, M. Grajcar, A. Korotkov, 
A. Shnirman, and A. Smirnov for a detailed discussion of 
the problem. We thank also A. Izmalkov, W. Krech, A. 
Maassen van den Brink, V. Shnyrkov, Th. Wagner, and 
A. Zagoskin, for fruitful discussions. 

The authors acknowledge the support from D-Wave 
Systems. Ya. G. acknowledges partial support by the 



(a) 




-0.010 -0.005 0.000 

*/*-1/2 













f=9 MHz 




f=12 MHz 


(b) 


f=18 MHz 







-0.010 -0.005 0.000 

*/* -1/2 



FIG. 1: The dependence of the tank phase on the flux when 
under the influence of high frequency irradiation. / = F/h is 
the power of irradiation in frequency units. 
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